August 2003 
hep-ph/0308254 



Implications of the DAMA and CRESST experiments 
for mirror matter-type dark matter 



R. Foot 1 

School of Physics, 
University of Melbourne, 
Victoria 3010 Australia 



Mirror atoms are expected to be a significant component of the galactic dark matter 
halo if mirror matter is identified with the non-baryonic dark matter in the Universe. 
Mirror matter can interact with ordinary matter via gravity and via the photon-mirror 
photon kinetic mixing interaction - causing mirror charged particles to couple to or- 
dinary photons with effective electric charge ee. This means that the nuclei of mir- 
ror atoms can elastically scatter off the nuclei of ordinary atoms, leading to nuclear 
recoils, which can be detected in existing dark matter experiments. We show that 
the dark matter experiments most sensitive to this type of dark matter candidate 
(via the nuclear recoil signature) are the DAMA/Nal and CRESST/Sapphire experi- 
ments. Furthermore, we show that the impressive annual modulation signal obtained 
by the DAMA/Nal experiment can be explained by mirror matter-type dark matter for 
|e| ~ 5 x 10~ 9 and is supported by DAMA's absolute rate measurement as well as the 
CRESST/Sapphire data. This value of |e| is consistent with the value obtained from 
various solar system anomalies including the Pioneer spacecraft anomaly, anomalous 
meteorite events and lack of small craters on the asteroid Eros. It is also consistent 
with standard BBN. 
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The DAMA/Nal experiment [1, 2] has been searching for dark matter and has ob- 
tained some very exciting positive results which merit serious consideration. While they 
have interpreted their data in terms of weakly interacting heavy particles an alternative 
interpretation will be suggested here. 

In the DAMA/Nal experiment the target consists of 100 kg of radiopure Nal. The 
aim of the experiment is to measure recoil energy of the Na, I atoms due to interactions 
of dark matter particles with their detector. Due to the Earth's motion around the 
sun, the rate should experience a small annual modulation: 



According to the DAMA analysis [2], they indeed find such a modulation over 7 annual 
cycles at more than 6a C.L. Their data fit gives T = (1.00±0.01) year and t = 144±22 
days, consistent with the expected values. [The expected value for to is 152 days (2 
June), where the Earth's velocity, ve, reaches a maximum with respect to the galaxy]. 
The strength of their signal is A = (0.019 ± 0.003) cpd/kg/keV. 

The DAMA collaboration have interpreted these impressive results as evidence for 
heavy weakly interacting dark matter particles. However, another possibility is that 
this experiment has observed the impacts of galactic mirror atoms, as will shortly be 
explained. 

Mirror matter is predicted to exist if nature exhibits an exact unbroken mirror 
symmetry [3] (for reviews and more complete set of references, see Ref.[4]). For each 
type of ordinary particle (electron, quark, photon etc) there is a mirror partner (mirror 
electron, mirror quark, mirror photon etc), of the same mass. The two sets of particles 
form parallel sectors each with gauge symmetry G (where G = SU (3)®SU (2) ® U(l) in 
the simplest case) so that the full gauge group is G®G. The unbroken mirror symmetry 
maps x — > —x as well as ordinary particles into mirror particles. Exact unbroken time 
reversal symmetry also exists, with standard CPT identified as the product of exact T 
and exact P[3]. 

Ordinary and mirror particles can interact with each other by gravity and via the 
photon-mirror photon kinetic mixing interaction: 



where F^ v (F' ) is the field strength tensor for electromagnet ism (mirror electromag- 
netism) 2 . Photon-mirror photon mixing causes mirror charged particles to couple to 
ordinary photons with a small effective electric charge, ee[3, 7, 8]. Interestingly, the exis- 
tence of photon-mirror photon kinetic mixing allows mirror matter to explain a number 
of puzzling observations, including the Pioneer spacecraft anomaly[9, 10], anomalous 

2 Given the constraints of gauge invariance, renomalizability and mirror symmetry it turns out[3] 
that the only allowed non-gravitational interactions connecting the ordinary particles with the mirror 
particles are via photon-mirror photon kinetic mixing and via a Higgs-mirror Higgs quartic interaction, 
C = \$ 4>cf)^ <p' . If neutrinos have mass, then ordinary - mirror neutrino oscillations may also occur[5, 
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meteorite events [11, 12] and the unexpectedly low number of small craters on the aster- 
oid 433 Eros[13, 14]. It turns out that these explanations and other constraints[15, 16] 
suggest that e is in the range 

1(T 9 £ |e| ~ 5 x 1(T 7 . (3) 

More generally, mirror matter is a rather obvious candidate for the non-baryonic 
dark matter in the Universe because: 

• It is well motivated from fundamental physics since it is required to exist if parity 
and time reversal symmetries are exact, unbroken symmetries of nature. 

• It is necessarily dark and stable. Mirror baryons have the same lifetime as ordi- 
nary baryons and couple to mirror photons instead of ordinary photons. 

• Mirror matter can provide a suitable framework for which to understand the large 
scale structure of the Universe [17]. 

• Recent observations from WMAP[18] and other experiments suggest that the 
cosmic abundance of non-baryonic dark matter is of the same order of magnitude 
as ordinary matter Qf, ~ Qdark- A result which can naturally occur if dark matter 
is identified with mirror matter [19]. 

If mirror matter is identified as the non-baryonic dark matter, then the dark matter 
halo will consist of compact objects such as mirror stars and planets, as well as a mirror 
gas and dust component. Evidence for mirror stars arises from MACHO observations 
[20, 21] (and to some extent from the puzzling 'isolated' planets[22]) while the existence 
of close-in extrasolar planets can also be viewed as mirror matter manifestations [23]. 
The amount of material in compact form is probably less than 50% (coming from the 
MACHO upper limit). Thus we expect a dark matter halo with a significant gas/dust 
component containing mirror H', He' + heavier mirror elements. Assuming a local halo 
dark matter energy density of 0.3 GeV/cm 3 , then the number densities of A' = H', He' 
and heavier elements is then given by 

, 0.3 GeV _ 3 

n A > = U'—j 77 cm ' (4) 

M A t 

where ^a> = Pa> / (0.3 GeV/cm 3 ) is the A' proportion (by mass) of the halo dark matter. 
As discussed above, a plausible value for Y^a'^a 1 is ~ 1/2. 

Arguments from early Universe cosmology (mirror BBN)[17] suggest that He' domi- 
nates over H', quite unlike the case with ordinary matter. Mirror elements heavier than 
H', He' will presumably come from nucleosynthesis within mirror stars, qualitatively 
similar to the ordinary matter case. In the ordinary matter case, the galactic relative 
(mass) abundance of elements heavier than if, He (collectively called 'metals' in the as- 
trophysics literature) is estimated[24] to be roughly Z g ~ 0.02 =^ iuetais/iHe ~ 0.10. 
These heavier elements are made up primarily (> 90%) of O, Ne, N, C which have 
Ma/Mp ~ 16 ± 4, Z = 8 ± 2. Thus, oxygen provides an excellent 'average' for ordi- 
nary elements heavier than helium - except perhaps for iron [which is about 10 times 



less abundant (by mass) than oxygen]. In the case of mirror element abundances, 
we would expect a qualitatively similar picture, i.e. O' (and elements with nearby 
atomic number) should dominate the energy density after H', He', with a possible small 
Fe' contribution. Thus we need only consider four mirror elements: H' , He', O', Fe' 
(where O' stands for Oxygen and nearby elements). Of course, quantitatively, the ratios 
£,o'/£,He',£,Fe'/£,o' are quite uncertain because of the different initial values for He' /H' 
(coming from mirror BBN) and other different initial conditions. Although the pro- 
portion of the various mirror elements in the halo (gas/dust ratio etc) is uncertain, the 
mass scale is not a free parameter: mirror hydrogen, H' , is predicted to have exactly 
the same mass as ordinary hydrogen, i.e. Mh> = 0.94 GeV, mirror helium, He', has 
mass M He , = 3.76 GeV etc 3 . 

In an experiment such as DAMA/Nal, the measured quantity is the recoil energy, 
E R , of a target atom. The minimum velocity of a mirror atom of mass M A > impacting 
on a target atom of mass Ma is related to E R via the kinematic relation: 



Interestingly, most of the existing dark matter experiments are not very sensitive to 
mirror matter-type dark matter because v min [Eq.(5)] turns out to be too high. This 
is because they either use target elements which are too heavy (i.e. large Ma) or have 
a E R threshold which is too high. For example, the CDMS experiment uses Ge as the 
target material and has a threshold of 10 keV[27]. This means that v min « 1600 km/s 
(for He'). Although there would be no cutoff velocity at the galactic escape velocity 
for He' due to He' self interactions, the number of He' with such high velocities would 
be negligible. The existing experiments with the greatest sensitivity to light mirror 
elements are the DAMA/NaI[l, 2] and the CRESST/sapphire experiments[28]. Both 
of these experiments will be examined in detail. 

When a mirror atom (of mass Ma>, atomic number Z') encounters ordinary matter 
(comprised of atoms with mass Ma, atomic number Z) Rutherford scattering can occur, 
with center of mass cross section: 4 



where v cm is the center of mass velocity of the impacting mirror atom and M red = 
M A M A '/(M A + M A ') is the reduced mass 5 . In Eq.(6), F x (qr x ) (X = A, A') are the 
form factors which take into account the finite size of the nuclei and mirror nuclei. 

3 It is possible to construct mirror matter models with broken mirror symmetry, in which case the 
masses of the mirror particles need not be the same as their ordinary counterparts [25], but these 
models tend to be more complicated and/or less well motivated than the simplest case of unbroken 
mirror symmetry [26]. 

4 Note that unless otherwise stated, we use natural units where U = c = 1. 

5 Due to the screening effects of the atomic electrons, the cross section is modified (and becomes 
suppressed) at small scattering angles [6 S ~ l/(MA>vr ) with r ~ 1CT 9 cm]. 
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(q = (2MaE r ) 1 ^ 2 is the momentum transfer and rx is the effective nuclear radius). A 
simple analytic expression for the form factor, which we adopt in our numerical work, 
is[29]: 

FxWx) = 3 x e w ' 1 (7) 

qr x 

with r x = 1.14X 1 / 3 fm, s = 0.9 fm. 

This cross section, Eq.(6), can be expressed in terms of the recoil energy of the 
ordinary atom, E R , and lab velocity, v (i.e. the velocity in Earth rest frame): 

^ X (8) 



dE R E 2 R v 2 
where 

A = Wa F 2 A {qr A )F 2 A ,{qr AI ) (9) 

Note the l/E R dependence. It arises because the dark matter particles interact electro- 
magnetically (i.e. via exchange of massless photons). This is quite unlike the standard 
WIMP case and therefore represents a major difference between mirror dark matter 
and standard WIMP dark matter. 
The interaction rate is 



dE R 
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= }^N T n A , / — \v\d a v 

x J dE R k 

= 2^ N T n A'~^2 / ,| dv (10) 
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where is the number of target atoms per kg of detector 6 and f(v,VE)/k is the 
velocity distribution of the mirror element, A', with v being the velocity relative to the 
Earth, and v E is the Earth velocity relative to the dark matter distribution. The lower 
velocity limit, v m i n (E R ), is obtained from Eq.(5), while the upper limit, v max = oo, 
because of A' self interactions (as we will explain in a moment). 
The velocity integral in Eq.(10), 

i(e r) . r (ii) 

Jv min (E R ) k\v\ 

is standard (as it occurs also in the usual WIMP interpretation) and can easily be eval- 
uated in terms of error functions assuming a Maxwellian dark matter distribution [29], 
f(v,v E )/k = (nv 2 )- 3 / 2 exp[-(v + v E ) 2 /v$\, 

I{E R ) = ^ [erf(x + y) - erf(x - y)\ (12) 



6 For detectors with more than one target element we must work out the event rate for each element 
separately and add them up to get the total event rate. 



where 

XE^,^. (13) 

For standard non-interacting WIMPs, vq is expected to be in the (90% C.L.) range[30], 

170 km/s ~ v ~ 270 km/s. (14) 

In the case of a halo composed of H', He', heavier mirror elements and dust particles, 
there are important differences due to mirror particle self interactions. For example, 
assuming a number density of n# e / ~ 0.08 cm" 3 [c.f. Eq.(4)] the mean distance between 
He' -He' collisions is 1/ (n He ,cr elastic ) ~ 0.03 light years (using a elastic ~3x 10" 16 cm 2 ). 
One effect of the self interactions is to locally thermally equilibrate the mirror parti- 
cles in the halo. The He' (and other mirror particles) should be well described by a 
Maxwellian velocity distribution with no cutoff velocity. [He' do not escape from the 
halo because of their self interactions]. A temperature, T, common to all the mirror 
particles in the halo can be defined, where T = M A 'Vq/2 (of course, T will depend on 
the spatial position). One effect of this is that vq should depend on M A > with 



v (A') = v (HeyM He ,/M A , (15) 

Thus, knowledge of Vq for He' will fix v o for the other elements. We will assume that 
the halo is dominated by He' (which is suggested by mirror BBN arguments [17]), with 
Vq = v (He') in the range, Eq.(14). 

The Earth motion around the sun produces an annual modulation in y: 

y ~ y + Aycosu(t - t ) (16) 

where y^ = {ve)/vq and Ay = Ave/vo (for He', y^ ~ 1.06, Ay « 0.07). The parameter 
t turns out to be June 2, and u = 2n/T (with T = 1 year). Expanding I(E R ) [Eq.12] 
into a Taylor series [making the y dependence explicit, i.e. I(Eji,y) = I(Er)]: 

I(E R ,y + Ay cos u(t-t ))= I (E R ,y )+ Ay cos u(t- t ) (^-] (17) 

{x-yo) 2 _|_ e -(x+yo) 2 _ 



y=yo 
and 
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The net effect is an interaction rate 

dE R ~ = dE R ~ + dE R ~ 



(19) 



where 
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Clearly, galactic mirror atom interactions will generate an annual modulation, A cos uj(t- 



to), in the event rate coming from the component. 

To compare these interaction rates with the experimental measurements, we must 
take into account the finite energy resolution and quenching factor. The quenching 

factor relates the detected energy (E R ) to the actual recoil energy (E R ), 

E R = q A E R (21) 

and for the DAMA experiment, q^ a , Qi have been measured to be approximately 
q Na ~ 0.30, qi ~ 0.09 [1]. The energy resolution can be accommodated by convolving 

the rate with a Gaussian, with a res ~ 0.16 E R (from figure 3 of Ref.[31]). 

The DAMA collaboration give their results in terms of the residual rate in the 
cumulative energy interval 2-6 keV, where they find that 

A exp = 0.019 ± 0.003 cpd/kg/keV (22) 

This number should be compared with the theoretical expectation: 

1 3 

1 r- " 

where 
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with Ej = 2.0 keV + AE * j (j = 0, 1, 2, ...) and AE = 1.0 keV. 

We have numerically studied A th . We find that the O' contribution to DAMA/Nal 
dominates over the He' (H') contribution provided that io'/iwe' ~ 7x 10~ 4 (£,o>/£,h> ~ 
4 x 10~ 8 ). The reason for this is of course clear: the actual threshold recoil energy is 
6.7 keV for A' — Na interactions 7 , which means that the v min [from Eq.(5)] is: 

v min {H' - Na) = 2830 km/s 

v m i n (He — Na) = 795 km/s 

Vmin(O f - Na) = 290 km/s 

Vmin(Fe' - Na) = 166 km/s (25) 



Numerically we find that mirror atom interactions with Na dominate over / for recoil energies 
above the 2 keV software threshold. This reason for this is clear: the threshold velocity is much lower 
for interactions with Na which is because a) Na is a much lighter clement than I [c.f. Eq.(5)] and b) 
the quenching factor for Na is 0.3 (c.f. with 0.09 for I) which means that the actual recoil threshold 
energy is 6.7 keV for Na and 22 keV for /. Also note that the corrections due to the form factor, which 
were taken into account using the simple analytic expression, Eq. (7) , are reasonably small (~ 5%) for 
A' = O' but larger (~ 30%) for A' = F e'. 



Because v m i n {H' — N a) ^> vq, any H' contribution to the DAMA/Nal signal is expected 
to be very tiny and we will neglect it. The quantity v m i n (He' — Na) is also quite high 
which suppresses the He' contribution relative to O' and Fe' . 

Interpreting the annual modulation signal in terms of O' and Fe', i.e. setting 
A th = A exp , we find numerically that: 



where the errors denote a 3 sigma allowed range (corresponding to 0.010 < A exp < 
0.028). The best fit region will also be affected by systematic uncertainties in the 
quenching factors, form factors and astrophysical uncertainties [e.g. uncertainties in 
Uo(-A')]. These uncertainties will increase the possible parameter range, however a 
detailed investigation of these effects we leave for the future. 

Because of the different masses of the two components, O' and Fe', their relative 
contributions can potentially be determined by the differential recoil energy spectrum, 
A> . In figure 1 we examine the representative possibilities a) DAMA signal is dom- 
inated by <y [i.e. £ ' = 0.10, f A # = for A' ^ 0'\, b) DAMA signal is dominated by 
Fe' [i.e. £p e / = 0.026, £a' = for A' ^ Fe'}. The case where A t h is made up of approx- 
imately equal contributions from both O' and Fe', corresponding to £o ; — ^Fe 1 = 0.05 
is also given. However, again we point out that a careful study of systematic uncer- 
tainties will be necessary before any definite conclusions can be made about the ratio 

Besides the annual modulation effect, DAMA/Nal has also measured the absolute 
event rate. This rate will contain the signal [dRo/dE R , Eq.(20), convolved with a 
Gaussian to incorporate the detector resolution] plus any background contribution. 
An interesting point is that the cross section, Eq.(8), rises sharply (oc l/E R ) at low 
E R , and this effect may show up in the data. [In any case, we should check that our 
absolute rate from the signal does not exceed the measured absolute rate]. In figure 
2 we plot the absolute rate with parameters fixed by the annual modulation signal, 
Eq.(26). Also plotted is the measured rate obtained from Ref. [32]. Interestingly, the 
data does indeed show a sharp rise at low E R which is compatible with the parameters 
suggested by the annual modulation effect. The shape of the measured rate at low E R 
is nicely fitted by both O' and Fe' dark matter, but the normalization prefers O' over 
Fe' dark matter. However possible small systematic uncertainties such as calibration 
errors may be present: a 0.1-0.2 keVee calibration error would be enough to allow Fe' 
to fit the data at low E R . 

Implicit in our analysis is that the mirror atoms can reach the DAMA detector 
from all directions, without getting stopped in the Earth. The stopping distance of a 
mirror atom, A' (of energy E' = ^M^'V 2 ) in ordinary matter (of atomic number density 
n = p/Ma) can easily be evaluated from: 




(26) 




where E^ ax can be obtained from Eq.(5) and E™ n = 1/(2toMa) (due to atomic 

-g^r J ~ 10]. Eq.(27) can be solved to give the energy of 
the mirror atom after travelling a distance x through ordinary matter: 

where L is the stopping distance: 

M\M A ,v\ 
~ 8irpe 2 a 2 Z 2 Z' 2 lO 

where t> j is the initial velocity of the mirror atom. The stopping distance in earth for 
He', O' and Fe' can easily be obtained from the above equation, giving: 

L(He') ~ 10 7 km for |e| = 4 x 10" 9 , ^ > v min (He') ~ 795 fcm/s, 
L(O') ~ 5 x 10 4 A;m for |e| = 4 x 10" 9 , ^ > v min (0') ~ 290 fcm/s. 
L(Fe') ~ 3 x 10 3 A;m for |e| = 4 x 10~ 9 , u< ~ 200 fcm/s. (30) 

Since L(He'), L(0') are much larger than the Earth's diameter, the retarding effect 
of the Earth is relatively small and no large diurnal effect is expected (in agreement 
with DAMA observations[33]). Mirror iron, may lead to a possibly large diurnal effect. 
However, dark matter detection experiments depend on |e| v 7 ?^" while the stopping 
distance in earth, depends just on |e|. The significant uncertainty in the size of £a' 
implies corresponding uncertainty in e and hence L(A'). It is therefore still possible 
for the DAMA/Nal signal to be dominated by the Fe' component, without leading to 
any significant diurnal effect. Note that experiments with a lower threshold (and hence 
lower value of v m i n ) will have a much greater sensitivity to the diurnal effect, so this 
effect may show up in future experiments. 

Let us now consider implications of this interpretation of the DAMA signal for 
other experiments. The CDMS/Ge experiment [27] has searched for nuclear recoils due 
to WIMP-Ge elastic scattering. This experiment has a threshold energy of 10 keV and 
the quenching factor is assumed (but not measured!) to be 1. This experiment finds 
just 4 events satisfying their cuts with 10 keV < E < 20 keV for their exposure of 
10.6 kg-day. However because the target consists of the relatively heavy element, Ge, 
and the threshold is relatively high, 10 keV, the sensitivity of the CDMS experiment 
to light mirror elements, is completely negligible. Assuming |e| Co' / 0. 10 = 4.8 x 10~ 9 
(as suggested from the DAMA/Nal experiment, if O' dominates the rate), we find 
numerically that the number of O' induced events (above the 10 keV CDMS threshold) 
is much less than 1 for their exposure of 10.6 kg-day. 



(28) 



If there happens to be a significant Fe' component, then this may potentially be con- 
strained by CDMS/Ge experiment. In the case where Fe' dominates the DAMA/Nal 
experiment, then |e|i/^_p e //0.026 = 4.8 x 10~ 9 . Numerically, we find that this implies 



26 events in the 10 keV < E < 20 keV range for CDMS, for their 10.6 kg-day exposure 
(c.f. just 4 detected events). The low rate obtained by CDMS experiment suggests 
that Fe' does not dominate over O'. However given possible experimental uncertain- 
ties, the case of Fe' dominance is probably not completely excluded. For example, the 
quenching factor may turn out to be somewhat less than I. For example, a value of 
0.6 would reduce the expected number of events from 26 down to 5 events which is 
consistent with the data. 

Clearly experiments with a lower threshold than DAMA/Nal might potentially 
provide more stringent constraints. The only experiment with a lower threshold than 
DAMA/Nal is the CRESST/Sapphire experiment [28]. That experiment uses 262 g 
sapphire crystals (Al 2 3 ) as the target medium with a low detection threshold of 
E R (threshold) = 0.6 keV. These features make CRESST/Sapphire particularly sensi- 
tive to low mass dark matter particles such as He',0' (and even Fe'). Unfortunately, 
the CRESST experiment does not have enough statistics to be sensitive to the annual 
modulation due to the Earth's motion around the sun, nevertheless the shape and 
normalization of the measured energy spectrum provide useful information 8 . We now 
study in detail the implications of mirror matter-type dark matter for this experiment. 

In this experiment the quenching factor is assumed to be approximately equal to 1 
(however, again this has not been specifically measured) [28] . As with the DAMA/Nal 
experiment, the recoil spectrum, Eq.(10), needs be convolved with a Gaussian curve 
(with a res ~ 0A2A7AE res 9 ) in order to take into account the finite energy resolution 
of the detector, 
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dE' R (31) 



dE R Jo dE' R y/2n<j, 
The CRESST collaboration present their results in terms of the quantity, 



C 3 



= ~~/TTp I j£dE R (32) 
AE JEj dE R 



where Ej = 0.6 keV + AE * j (j = 0, 1, 2, ...) and AE = 0.2 keV. We have numerically 
studied Cj for various cases. In figure 3 we plot the expected value for Cj for the 
best fit values of |e| v / Ca' assuming the DAMA/Nal rate is dominated by a) O', b) Fe' 
and c) 50-50 O', Fe' mixture. [Recall, these are the same three cases which were fitted 

8 Because of the low threshold, the CRESST experiment might be sensitive to the diurnal effect 
and this could even show up in the existing data. 

9 Note that a res = j== ^ 0.4247A£' res implies a FWHM of AE res for the Gaussian curve, which 
is the CRESST prescription [28]. In Ref.[28], two values of AE res are discussed, AE res w 0.2 keV 
(from an internal calibration source) and AE res w 0.5 keV (from possible contamination with 55 Fe). 
In our numerical work we have used the former value (unless otherwise stated). Using AE res = 0.5 
keV would lead to |e|V6i' values smaller by about 20%. 



to the DAMA/Nal annual modulation signal and were plotted in figure 1]. While the 
shape of the CRESST/Sapphire data (obtained from figure 10 of Ref. [28]) is reasonably 
consistent with the expected shape from A' interactions, the normalizable is roughly a 
factor of 2 too high. This may be due to systematic uncertainties which we illustrate 
in figure 4. In this figure, the CRESST quenching factor is taken to be 0.7 instead of 
the assumed value of 1.0. (Similar results occur if there happens to be a small energy 
calibration uncertainty of ~ 0.2 keV). Figure 4 clearly demonstrates the rather nice 
fit of O', Fe' dark matter to the shape and normalization of the CRESST data (after 
allowing for reasonable systematic uncertainties). 

Given the rather nice fit of the shape and normalization of the CRESST data 
(within reasonable systematic uncertainties) to the expectations of Fe', O' dark matter 
from the DAMA/Nal fit, it is clearly very tempting to suppose that the CRESST data 
may be mostly signal with very little background component. On the other hand, 
the CRESST collaboration [28] have argued that their data is most likely background 
because of the rate of coincidence events. This argument required the background to 
be due to single particle interactions and isotropic which it may not be. 

Finally, note that the CRESST/Sapphire experiment is much more sensitive to 
H' , He' than the DAMA/Nal experiment. Assuming a pure He' halo, i.e. fa = for 
A' ^ He', we find that the CRESST data suggest: 



In this pure He' halo limit, the (CRESST) value for |e| v^£hv> above, is not consistent 
with the value from DAMA/Nal [Eq.(26)]. Thus He' cannot dominate the rate for 
DAMA or CRESST. This suggests that: fa ~ 0.2^ and/or £ Fe , ~ 0.04^ ■ Clearly 
this constraint is significant, but nevertheless, still allows He' to be the dominate halo 
dark matter component 10 . 

Assuming that DAMA and CRESST have detected galactic mirror matter-type 
dark matter, then this suggests an |e| value of around 10 -8 — 10~ 9 . Previous work 
(see Ref. [13] and references there-in) looking at various solar system implications of 
mirror matter has identified a similar but somewhat larger range for e, Eq.(3). This 
information is summarized in figure 5. Also shown is the experimental bound[16, 15], 
lei ~ 5xl0" 7 coming from recent orthopositronium lifetime measurements [34] and also 
the limit suggested from BBN[35]. 

Let us also mention that if |e| ~ 5 x 10~ 9 , there will be interesting terrestrial effects 
of mirror matter. Fragments (of size R) of impacting mirror matter space bodies can 
remain on/near the Earth's surface provided that [36] 



10 Note that the CRESST/Sapphire experiment does not put significant limits on the H' component 
of the dark matter halo. Numerically, we find that He' dominates over H' provided that £jy <~ 15£hv 
which is not a very stringent condition. 




(33) 




(34) 



Such fragments can potentially be detected and extracted with a centrifuge [36]. If 
mirror matter fragments become completely embedded within ordinary matter (which 
is necessarily the case for e < 0) then the fragments will thermally equilibrate with the 
ordinary matter environment. The observational effect of this is to cool the surrounding 
ordinary matter, as heat is transferred to the mirror body and radiated away into mirror 
photons[37]. Finally, even tiny solar system mirror dust particles can lead to observable 
effects. These particles impact with the Earth with velocity in the range 11 km/s ~ 
v ~ 70 km/s and can be detected in suitably designed surface experiments [38] such as 
the St. Petersburg experiment [39]. 

In conclusion, we have pointed out that the DAMA/Nal, CRESST/Sapphire and 
other dark matter experiments are sensitive to mirror matter-type dark matter. Fur- 
thermore, the annual modulation signal obtained by the DAMA/Nal experiment can 
be explained by mirror matter-type dark matter for |e| ~ 5 x 10~ 9 . This explanation 
of the DAMA signal is supported by DAMA's absolute rate measurement as well as by 
the size and shape of the CRESST data. Furthermore this explanation is not in con- 
flict with CDMS or any of the other dark matter experiments because of their higher 
thresholds. The e value suggested by the DAMA/Nal experiment is consistent with the 
value obtained from various solar system anomalies including the Pioneer spacecraft 
anomaly, anomalous meteorite events and lack of small craters on the asteroid Eros. 
It is also consistent with standard BBN. 
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Figure 2: The absolute event rate for the DAMA/Nal experiment for O' , Fe' dark matter (for 
vq = 230 km/s). The parameters are given by the fit to the DAMA/Nal annual modulation 
effect, where we take the same three representative cases as figure 1: a) Hv^O'/O-lO = 
4.8 x 10~ 9 , U> = for A' / O' (short-dashed line) b) |e|V&v/0.026 = 4 -8 x 10~ 9 , £ A , = 
for A' / Fe' (long-dashed line) and c) |e| = 4.8 x 10" 9 with = 4£ Fe / = 0.05 (dotted line). 
Also shown is the DAMA/Nal data obtained from Ref.[32]. 
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Figure 5: Favoured range of e from various experiments/observations. 



